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We demonstrate that a completely integrable classical mechanical model, namely the lattice 
Landau-Lifshitz classical spin chain, supports diffusive spin transport with a finite diffusion con- 
stant in the easy- axis regime, while in the easy-plane regime it displays ballistic transport in the 
absence of any known relevant local or quasi-local constant of motion in the symmetry sector of the 
spin current. This surprising finding should open the way towards analytical computation of diffu- 
sion constants for integrable interacting systems and hints on existence of new quasi-local classical 
conservation laws beyond the standard soliton theory. 



Introduction.- Derivation of irreversible macroscopic 
transport (e.g. Fourier/Ohm/Fick's) laws from re- 
versible, deterministic, microscopic equations of motion 
is one of the central questions of statistical physics which 
remains largely unsolved even today. It has been believed 
[IHB] that chaotic dynamics in classical systems, or more 
generally strong non-integrability in either quantum or 
classical systems, are necessary conditions for diffusive 
transport. Recently, few examples of spin-diffusion at 
high temperature in completely integrable but strongly 
interacting quantum spin/particle chains have appeared 
[7HTT1. suggesting that complete integrability might not 
exclude the possibility of macroscopically diffusive dy- 
namics. It has remained unclear, however, whether quan- 
tum nature of the corresponding many-body dynamics 
supporting macroscopic entanglement is a necessary con- 
dition. Here we show that even quantum correlations are 
not necessary. By performing extensive numerical simu- 
lations in a family of integrable classical spin chains with 
local interactions - the lattice Landau-Lifshitz model [T2] 
- we show that spin transport at finite temperature is dif- 
fusive in the easy axis regime, while it becomes ballistic 
in the easy-plane regime and anomalous at the isotropic 
point. This opens up the possibility for analytic compu- 
tations of diffusion constants in interacting many-body 
systems. 

Liouville integrability [13 is the central concept in the 
analytic theory of classical mechanics. A Hamiltonian, 
i.e. conservative system in classical mechanics, is inte- 
grable if it possesses the same number of independent 
conserved quantities as the number of degrees of free- 
dom, call it n. In other words, its motion can be re- 
duced to quasi-periodic winding around n— dimensional 
torus embedded in 2n dimensional phase space [13 . In- 
tegrable dynamics is thus regular and manifestly free of 
sensitive dependence on initial conditions. Nevertheless, 
integrable systems, though being sparse in nature, repre- 
sent one of the key topics in mathematical physics as they 
gave birth to the celebrated soliton theory [12 explaining 



a variety of observable phenomena, ranging, to name just 
a few, from localized light in non-linear optics, waves on 
shallow water and tsunami waves, to elementary parti- 
cles and localized excitations in condensed matter at low 
temperatures. 

The solitons, indestructible localized packages of en- 
ergy which propagate through the system and scatter 
from each other like elastic hard balls, have been be- 
lieved to be the reason why integrable extended systems 
behave as ideal - ballistic conductors of heat, particles, 
electric charge, magnetization, etc. (TJ [2]. Being par- 
ticularly interested in the one-dimensional lattice sys- 
tems, where n particles are arranged along a line or a 
ring such that only nearest neighbors can interact rep- 
resenting the simplest model of crystalline solids, one 
finds that existence of non-trivial conservation laws (be- 
sides the transported quantity, e.g. energy, particle num- 
ber, electric charge, magnetization) generically implies 
the ballistic (non- diffusive) transport [H]. This state- 
ment, which builds on an old idea of Mazur [TS] but has 
only recently been formally proven [16 , essentially states 
the following. Whenever there exist a quantity / which 
is conserved in time L{t) = / and independent of the 
transported quantity itself, such that (/J) ^ 0, where 
Jif) — Xl^=ii(^5^) is ^^^ current with j{x^t) being the 
current density at time t and at site x in the lattice and 
(•) denotes the thermodynamic average (for fixed, spec- 
ified values of temperature, electro-chemical potential, 
magnetization, etc.), then the transport is ballistic and 
the corresponding Kubo conductivity hi diverges. Con- 
ductivity is related to a diffusion constant D^ via the 
generalized Einstein relation K: = D/T where T is the ab- 
solute temperature, and the latter can be within Green- 
Kubo linear response theory expressed 



D = lim lim 



>oo n^^oo 



n .r 



C{r, t)dt, 



(1) 



in terms of the integrated spatio-temporal current- 
current correlation function C(r, t) = (j(x, 0)j{x + r, t)). 



Lattice site x is arbitrary for translationally invariant 
systems which exhaust major examples of integrable sys- 
tems. Note that in statistical mechanics the thermody- 
namic limit of size n ^ oo has to be taken always before 
the time r ^ oo limit. It is clear that the above expres- 
sion for the diffusion constant can be given in terms of 
time-correlations only, namely D = \im.t^oo D{t), intro- 
ducing a time-dependent diffusion constant as 



D{t) 



C{t)dt 



(2) 



where we can write the total current time autocorrela- 
tion as C{t) = -(J(O)J(t)). The existence of solitons 
and nontrivial conserved quantities in integrable systems 
implies non- vanishing tails of the time correlations [HI [16] 
C'oo = C(t ^ oo) ^ 0, in turn implying linear divergence 
of the time-dependent diffusion constant D -^ Coot -^ oo 
which is a signature of ballistic transport. Up to date, 
all studies of transport in classical integrable particle 
chains have persistently showed ballistic transport (see 
e.g. Ref. pT and references therein). On the other hand, 
for diffusive transport {D < oo), we have evidence that 
not even microscopic chaos [6 is neccessary but a weaker 
property of dynamical mixing [18 is sufficient p!7|[T9l [2Q]. 
Furthermore, some recent numerical studies of quantum 
spin 1/2 chains with anisotropic Heisenberg interaction 
{XXZ chains) indicated [1 H [HI [22] that the high 
temperature spin transport is diffusive in the easy-axis 



regime, despite the fact that the XXZ chains are quan- 
tum integrable by the algebraic Bethe ansatz which is the 
quantum version of the soliton theory. These results have 
been further corroborated with evidence of particle and 
spin diffusion in another Bethe Ansatz integrable model 
namely the one-dimensional Hubbard model [11]. Nev- 
ertheless, as the completeness of Bethe ansatz solutions 
has not been proven in these models and as quantum dy- 
namics supports a high degree of complexity as opposed 
to classical dynamics due to exponentially large (in n) 
Hilbert space dimension and possibility of macroscopic 
entanglement [23], it has remained unclear whether the 
key to the observed diffusion really lies in the integrabil- 
ity structures of the Heisenberg and Hubbard models. 

Model and methods.- However, in this Letter we show 
a convincing numerical evidence for spin diffusion at high 
temperature in a fully classical completely integrable 
model. We consider a lattice Landau-Lifshitz (LLL) 
model [12] for a chain or a ring of n classical spins de- 
scribed by angular-momentum vectors 6^, x = 1, ... n of 
equal fixed length |6^| = R and the Hamiltonian function 
of the form 



H — 2_^ h{Sx, Sx-\-i), 



(3) 



where the nearest-neighbor spin interaction density 
h{S^ S') is given by 



h{S, S') = log |cosh(p^3) cosh{pS'^) + coth^ {pR) sinhipSs) smh{pS'^) + sinh-^ {pR)F{Ss)F{S'^){SiS[ + ^s^^)] (4) 



and F{S) = ^J{smh^{pR) - sinh^ {pS)) / {R^ - S'^) (as defined by Faddeev and Takhtajan, Ref.[l2], Chapter III.5). p 
is the model's parameter, which can be real or purely imaginary. So we shall reparametrize it with an alternative 
real parameter 5 = p^ ^ which may - in analogy with the closely related XXZ spin chain [24 - be called anisotropy 
parameter. The cases with 6 > correspond to easy- axis spin-spin interaction, those with S < to easy-plane 
interaction, whereas S = designates the case of isotropic interaction. 

As the third component Ms of the total magnetization M = XI ^=i '^x is a constant of motion, one finds that the 
equation of motion for local magnetization dSx^s/dt = jx — jx-i has a form of a local continuity equation from which 
we read out the expression for the spin current density jx = j{Sx^ Sx-\-i) in LLL model, namely 



j{Sj') = smh-\pR){S2S[ - SiS'^)F{Ss)F{S^s)^M-HSj') 



(5) 
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The Hamiltonian and the current density in the easy- 
plane regime ^ < can be obtained from (4|5 ) by analytic 
continuation, i.e., by replacing cosh(p») by cos(5») and 
sinh(p») by sin(5») where s = ^/^. In the isotropic case 
^ = one performs the limit p — ^ in (ip) and obtains 
explicitly 



h{S J') = log [1 



S-S' 



.j{sj') 



020 1 — S1S2 

R^ + S-S' 



(6) 



Time development of each spin is given by solving Hamil- 



ton's equations dSx/dt = ^^ x Sx^ where x denotes the 

cross product. These equations of motion for Sx{t) are 
then solved using a variable step integrator of Matlab 
with required relative accuracy of 10~^ for all trajecto- 
ries. Besides checking the accuracy of trivial conserva- 
tion laws, such as H and M3, we have also checked that 
numerically determined Lyapunov exponents [13 vanish 
asymptotically in all three regimes {S >,=,< 0) as re- 
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FIG. 1. (Color online) Current autocorrelation function C{t) 
in log-log scale for easy-plane regime (top curves, orange: 
n = 160, black: n = 2560), isotropic regime (middle curves, 
yellow: n = 2560, blue: n = 5120) and easy- axis regime 
(bottom curves, violet: n = 2560, green: n = 5120). Shaded 
regions denote the estimated statistical error for ensemble av- 
erages over N ^ 10^ initial conditions. Dashed lines denote 
asymptotic behavior for large time in the easy-plane regime 
(dark-blue) and isotropic regime (light-blue). 
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FIG. 2. (Color online) Time dependent diffusion constant 
D(t) in isotropic regime (left) and easy- axis regime (right) 
comparing different system sizes n = 640 (black), n = 2560 
(orange) and n = 5120 (blue). The shaded region denotes es- 
timated statistical error and the dashed hues show estimated 
asymptotic behavior. 



quired for a completely integrable system. 

We chose N initial conditions generated by the 
Metropolis algorithm [25] yielding a thermal Gibbs en- 
semble with given temperature T and vanishing Ms = 
component (along the symmetry axis of the spin-spin in- 
teraction) of the magnetization. We note that ensembles 
with non-vanishing or non-fixed Ms would not yield the 
correct Kubo formula [4] for the spin-diffusion constant 
in the absence of external magnetic field. The LLL model 



is invariant under a canonical transformation 

i.e. a TT-rotation around the first cartesian axis, namely 
TZH = H. As the LLL r— matrix [12] is invariant under 
IZ so are also all the derived local conserved quantities 
IZIk = Ik- However, the spin current is odd under the 
symmetry transformation, 7ZJ = —J, so all (ballistic) 
terms in the Mazur bound [14] have to vanish, {Ik J) = 
{TZIklZJ) = —{Ik J), allowing for the possibility of non- 
ballistic transport. 

Results.- Solving Hamilton's equations numerically we 
compute the spatio-temporal current-current autocorre- 
lation function C(r, t), as well as the temporal autocorre- 
lation C{t) and obtain the time-dependent spin diffusion 
constant D{t) ^. We chose several different chain sizes 
n (up to n = 5120) and ensemble size N large enough 
(typically N ~ 10^) so that finite size effects, and the 
statistical error scaling as ~ I/a/tl/V, appear negligible. 
Representative values of the anisotropy parameters are 
chosen as (5 = 1 (easy- axis), S = —1 (easy-plane), and 
5 = (isotropic regime) while fixing R = 1. In all 
cases the temperature of the initial Gibbs ensemble is 
T = 4. In Fig. [l] we plot spin current autocorrelation 
C{t) in all three regimes and find a finite saturation value 
C'oo 7^ implying ballistic transport in the easy-plane 
regime, whereas, in the isotropic and easy-axis regimes 
the tails of the current autocorrelation function are van- 
ishing. To elaborate on the tails of the current autocor- 
relation function we plot in Fig. [2] the time dependent 
diffusion constant D{t)^ which in the isotropic regime 
shows power low behavior D{t) ex t" with a ~ 0.35 
and in the easy-axis regime saturates at a finite value 
D{oo) ^ 0.38, providing a firm evidence of anomalous 
spin-transport in the isotropic regime and diffusive spin- 
transport in the easy-axis regime. These results can be 
even better illustrated by portraying spatio-temporal cor- 
relations (7(r, t) in Fig. [s] In the easy-plane regime we 
find a clear causality-cone structure with non-decaying 
tails whereas in the easy-axis case the cone is curved in- 
wards and the tails are strongly damped in time result- 
ing in suppression of all ballistic contributions for long 
times. In bottom panel of Fig. [3] we present data for 
much smaller ensemble of initial conditions and smaller 
system size showing "scars", i.e. ballistic soUtons which 
may emerge in initial conditions with localized large ther- 
mal fluctuations and which can travel much faster than 
the correlation velocity, but contributions of which (due 
to the existence of IZ symmetry) vanish in the limit of 
infinitely large ensemble N ^ oo. Therefore, one would 
recognize the fact that we are dealing with a completely 
integrable system (in particular in the easy axis regime) 
only by looking at a finite N data. 

Conclusion.- The surprising possibility of having a 
regime with normal, diffusive transport in completely in- 
tegrable classical system (as found in the easy- axis regime 



of LLL model), which is associated with the existence 
of parity-type symmetry of the classical r-matrix whose 
operation changes the sign of the corresponding trans- 
porting current, opens up the immediate question if the 
classical r-matrix theory fT2[ can be updated to allow for 
analytical calculations of diffusive transport coefficients. 
Equally interesting is the problem of explaining the bal- 



listic regime in such a situation (say in the easy-plane 
regime of LLL model), since the Mazur bound identically 
vanishes. As LLL model can be considered as a classical 
limit of XX Z spin chain (see also [21]), we conjecture 
the existence of a classical analog of the construction of 
quasi-local conserved quantities with negative parity [27] 
going beyond the standard soliton theory [12] . 
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FIG. 3. (Color online) Modulus of spatiotemporal spin current autocorrelation function |C(r, t)| shown in log-scale with color 
scale ranging from 10~^"^ to 10~^ indicated in the bottom-right. In the upper panels we show data averaged over ensembles of 
N ^ 10^ initial conditions in easy-axis (left; n — 5120), isotropic (center; n — 5120 ) and easy-plane (right; n — 2560) regimes. 
The bottom plot shows easy-axis regime again with a smaller system/ensemble size n — 160, N — 600 where scars of solitons 
emerging from local thermal fluctuations are still clearly visible. 
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